Avalanches and clusters in planar crack front propagation 
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We study avalanches in a model for a planar crack propagating in a disordered medium. Due 
to long-range interactions, avalanches are formed by a set of spatially disconnected local clusters, 
the sizes of which are distributed according to a power law with an exponent r a = 1.5. We derive 
a scaling relation r a = 2r — 1 between the local cluster exponent r a and the global avalanche 
exponent r. For length scales longer than a cross-over length proportional to the Larkin length, the 
aspect ratio of the local clusters scales with the roughness exponent of the line model. Our analysis 
provides an explanation for experimental results on planar crack avalanches in Plexiglas plates, but 
the results are applicable also to other systems with long-range interactions. 

PACS numbers: 62.20.mt, 45.70.Ht, 64.60. av 



I. INTRODUCTION 

Driven elastic manifolds in disordered media have been 
used to model a number of physical systems ranging from 
crack propagation in solids [l|-|8jl to dynamics of magnetic 
domain walls in ferromagnets |9| and vortices in type-II 
superconductors [l(J, and to charge density waves [TTj | . 
Such systems exhibit non-equilibrium phase transitions 
as an external driving force / is varied, so that below 
the depinning threshold f c the system is pinned by the 
disorder, while for f > f c the system moves at a finite 
average velocity. For / close to but above f c , such motion 
typically occurs as a sequence of avalanches with a broad 
distribution of sizes. A large class of such systems, in- 
cluding advancing crack fronts in solids 0, , contact 
lines in wetting [12| , magnetic domain walls with dipolar 
interactions [ij, and plastically deforming crystals [l3j . 
is characterized by a long-range interaction kernel. The 
peculiarity of these systems is that due to the long-range 
interactions governing the avalanche dynamics, the area 
swept over by an avalanche is not necessarily simply con- 
nected in space as it would be the case if the interactions 
were only short-ranged. Instead, global avalanches are 
formed by a set of spatially disconnected local clusters. 

The propagation of a planar crack confined in a weak 
plane represents an ideal system to realize experimen- 
tally the depinning transition predicted by the crack line 
model Yet, the first experimental results for the 

crack front roughness 0, H[ and for the avalanche size 
distribution [6j disagreed with the theory. The rough- 
ness exponent of the crack front was found in the range 
£ = 0.55 — 0.6 0, Q, instead of the theoretical value 
C = 0.39 [3, while the avalanche sizes defined by the 
waiting time matrix measuring only spatially con- 
nected parts of the avalanches, were found experimen- 
tally to be power law distributed with an exponent much 
higher than the value predicted by the long-range line 
model for avalanches extracted by considering the global 
velocity signal of the crack front. These earlier claims 
have been revised recently: by testing a larger range of 



length scales it was shown that the roughness exponent 
would cross over to the theoretical value at sufficiently 
large scales [l5|,[l6j]. Furthermore, in Ref. Q it was shown 
that the crack line model can reproduce the avalanche 
statistics of the experiment, if also the numerical data 
are analyzed using the waiting time matrix introduced 
in while a smaller exponent value is obtained when 
considering avalanches extracted from the global velocity 
signal of the crack front model. Still, a complete under- 
standing of the origin of the different scaling exponents 
is lacking. 

Here we clarify these issues by studying the statistics 
of global avalanches and their spatially connected parts, 
the local avalanches (or clusters), for a long range elastic 
string moving in a disordered medium under quasi-static 
external driving. Both the avalanche and cluster size 
distributions are found to exhibit scaling, but with dif- 
ferent power law exponents. In particular, by studying 
the model with a larger range of length scales, we obtain 
a more accurate description of the scaling behavior, and 
therefore a better estimate of the related exponents than 
in previous numerical studies 0] ■ These numerical results 
are compared with data from experiments on Plexiglas 
plates [6J, and excellent agreement is found. Moreover, 
we propose a scaling relation between the power law ex- 
ponents of the avalanche and cluster size distributions by 
considering the avalanche break-up process, supported by 
our numerical results. We finally study the morphology 
of the local clusters, and find that their aspect ratio scales 
with the roughness exponent of the line model above a 
cross-over length proportional to the Larkin length of the 
crack front. Previous studies [7] using a continuous time 
version of the crack line model have presumably been 
sampling the short length scale regime. Thus, our re- 
sults explain and clarify a number of open issues related 
in particular to planar crack front propagation, such as 
the different scaling of global avalanches and local clus- 
ters, as well as the connection between roughness and 
avalanche morphology. However, the implications of our 
study extend beyond crack front propagation, as the re- 
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suits can be immediately extended to a number of other 
systems described by the same model, such as contact line 
dynamics in wetting JTa , and low angle grain boundaries 
in crystalline solids [17] . Furthermore, a number of other 
driven systems in which the avalanche dynamics is gov- 
erned by long range interactions, ranging from domain 
walls in ferromagnets with dipolar interactions [t| to plas- 
tically deforming crystalline solids [l3[, are expected to 
obey similar scaling laws. The paper is organized as fol- 
lows: In the next Section, the simulation model and the 
methods used to extract the local clusters from the sim- 
ulations and the experiments are presented. These are 
then followed by the results of the numerical simulations 
along with a direct comparison with experimental results 
from Plexiglas plates in Section III. Section IV finishes 
the paper with conclusions and discussion. 



II. MODEL AND METHODS 
A. Quasi-static model of crack front propagation 

The model of the propagating planar crack front, in- 
troduced in Ref. @ , is represented by a vector of integer 
heights hi, i = 1,...,L, with L the system size. Crack 
propagation is driven by the local stress intensity factor 
(SIF) Ki, which represents the asymptotic prefactor of 
the 1 / y/r divergence of the stress field near the crack tip. 
To model the different contributions of the local SIF Ki 
acting on a front element i, it is taken to be of the form 
Ki = Kf astlc + Kl a h ™ dom + K ext , where 
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FIG. 1: a) An example of the avalanche map A(x,y), for 
k = 0.05 and To = 1, with the area swept over by differ- 
ent avalanches denoted by different shades of gray (time is 
running from black to white), b) The same data as in a) 
represented by the velocity matrix V{x, y) = 1/W(x,y) with 
a zero threshold. The crack front is moving to the positive y 
direction. Black corresponds to zero velocity, while regions of 
finite velocity are shown in white, c) An example of a struc- 
ture of a single avalanche as identified by the AM method, 
d) The same avalanche as in c) extracted by using the WTM 
method. 
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represents the first order variation of the stress intensity 
factor due to a first-order perturbation of the front po- 
sition, b is the front segment spacing and Tq tunes the 
strength of the elastic interactions [l[ , K r i a £ dom is a time- 
independent Gaussian random variable of zero mean and 
unit variance, with 
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representing random toughness of the material, and K ext 
is the contribution of the external load. Notice that the 
long range kernel in Eq. ([I]) is identical in many other 
problems, such as contact lines [12[ and low angle grain 
boundaries [ijj ■ Periodic boundary conditions are im- 
posed. The dynamics is defined in discrete time t by 
setting 



v i (t) = h i (t + l)-h i (t) = e(K i ), 



(3) 



where Vi is the local velocity of the front element i, and 
9 is the Heaviside step function. During a single time 
step, the front elements with Vi > are advanced by a 
unit step, hi(t + 1) = hi(t) + 1, and new random forces 



are generated for those elements. The local forces are 
then computed again for each element, and the process 
is repeated until Vi = for all i and the avalanche stops. 
The size s of the avalanche is the total number of elemen- 
tary moves during such an avalanche. Then the external 
stress is increased so that exactly one of the elements be- 
comes unstable, and a new avalanche is initiated. While 
such a discretization of the dynamics neglects the fact 
that the local velocity of the crack front should be pro- 
portional to the local SIF, this is a standard technique 
employed in models of the depinning transition, going 
back to the work of Leschhorn [18j , and is known to have 
no influence on the scaling behavior we study here. More- 
over, this procedure is essential to be able to drive the 
system quasi-statically. Such quasi-static driving has the 
advantage over continuous time models with a finite driv- 
ing velocity [7] that avalanches and clusters can be de- 
fined without ambiguity also for a zero threshold. As the 
crack front advances, the applied SIF K ext decreases at 
a rate proportional to the instantaneous average velocity 
v(t) = l/LJ2i=i v i(t) °f the front, with a proportional- 
ity constant k. Thus, the cut-off of the avalanche size 
distributions is expected to scale with k Q. 
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Methods to extract the local clusters 



To extract the clusters from an avalanche, we use two 
methods. The first is the avalanche map (AM), based on 
a two dimensional array A(x, y) such that a unique value 
is assigned to all points (x, y) over which an avalanche has 
swept. We then define the local cluster size a AM as the 
area of spatially connected regions with the same value of 
A(x, y) (see Fig. Q]a and c). The second approach is the 
recently proposed waiting time matrix (WTM) W(x,y), 
defined by the time the interface has spent within a pixel 
corresponding to the location (x,y) [1, 0- The local 
velocity matrix is then given by V(x,y) = 1/W(x,y). 
Clusters are defined as spatially connected regions of 
area a WTM within which the local velocity exceeds some 
threshold value vth (see Fig. [T]b and d). Here we study 
the model in the quasi-static limit, which allows us to 
use effectively a zero threshold, by setting the waiting 
time of the pinned configurations separating avalanches 
to be much longer than the maximum avalanche dura- 
tion. When using the WTM approach to record the mo- 
tion of the crack front during an avalanche, only n — 1 
steps are recorded when a line element moves n steps, as 
both the initial and final pinned configurations are clas- 
sified as immobile. Thus, the avalanche in Fig. [T]d) is a 
bit thinner in the y direction as compared to the one in 
Fig. He), and parts of the avalanche involving motion of 
a single step forward only are lost completely. However, 
our results suggest that this docs not affect the statistical 
properties of the clusters. In the case of the discrete line 
model, another issue is that for soft lines (Tq < 0.75), 
in which the height difference between neighboring seg- 
ments of the line can be two or more pixels of the WTM, 
two consecutive avalanches can be recorded as one. In the 
following we use WTM only for large enough Tq . Notice 
that the AM does not suffer from such problems. 



III. RESULTS 

A. Avalanche and cluster size distributions 

The probability distributions of avalanches and cluster 
sizes for L = 2048, T = 1 and various values of k are 
shown in Figs. [5] and [3J respectively. The sizes s of 
the avalanches (Fig. [2] open symbols) are distributed 
according to a power law with a fc-dependent cut-off, 



P(s) = s- r f k (s/k- 1 ^), 



(4) 



with fk(x) a scaling function describing the shape of the 
cut-off, r = 1.25 ± 0.05 and l/a k = 0.725 ± 0.08. The 
value of r is in good agreement with results from simu- 
lations and unpublished experiments |19j . However, 
the scaling of the cut-off was not considered in Ref. Q. 
Notice that the scaling function f k {x) displays a bump 
around the cut-off of the distribution [2(|, and thus the 
true power law exponent can be seen directly only for 
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FIG. 2: (color online) The main figure shows the size dis- 
tributions of the avalanches for To = 1 and various values 
of k (open symbols), as well as the size distributions of the 
avalanches for a fixed k — 0.0125 and various values of To 
ranging from To = 0.4 to To = 1.0 (filled symbols, with 
the latter set of distributions displaced vertically for clar- 
ity). The solid line is a guide to the eye and corresponds 
to t = 1.25. The inset shows a data collapse of the distribu- 
tions with To = 1 and various values of k, with r = 1.25 and 
l/cr fe = 0.725. 



s well below the cut-off scale. The exponent values we 
quote are obtained by optimizing the data collapse of 
the distributions. To demonstrate the robustness of our 
results with respect to changing the values of various pa- 
rameters of the model, we also show the avalanche size 
distributions for a fixed k but varying Tq in Fig. [2] (filled 
symbols): the distributions turn out to be virtually inde- 
pendent of Tq in the range 0.4 < Tq < 1.0 considered in 
Fig. [2] Analogously, the cluster sizes or areas a AM (top 
panel of Fig. |3j and a w ™ (bottom panel of Fig. |3j 
extracted by the AM and WTM methods, respectively, 
are both observed to scale according to 



P{a) = a- T «g k (a/k- 1 ^x), 



(5) 



but with an exponent r a significantly different from r, 



r AM 



1.52 ±0.05 and r, 



WTM 



1.53 ±0.05. Thus, 



the two methods to define the clusters give the same r a - 
exponent within errorbars. In this case the scaling func- 
tion g k (x) does not exhibit a bump, and the scaling of 
the distributions extends all the way to the cut-off scale. 
Also for the cluster size distributions similar results are 
obtained for different values of Tq (not shown). The 
observed exponent value is somewhat lower than that 
quoted in Refs. [1, 0] for the WTM approach, but our 
analysis below indicates that also the data presented in 
d, 0] is perfectly consistent with t^™ rj 1.5. This is 
because in Refs. 0,0] also the region of the cut-offs of the 
distributions was included in the fits, thus overestimat- 
ing the true exponent value. Notice also that the cut-off 
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FIG. 3: (color online) Top: The size distributions of the AM 
clusters for To = 1 and various values of k. The solid line 
is a guide to the eye and corresponds to r a = 1.52. The 
inset shows a data collapse with r a = 1.52 and l/o"& = 0.7. 
Bottom: The size distributions of the WTM clusters for To = 
1 and various values of k. The solid line is a guide to the eye 
and corresponds to r a = 1.53. The inset shows a data collapse 
with r a = 1.53 and l/<7fc = 0.8. 



FIG. 4: (color online) The scaled size distributions of the 
WTM clusters for various threshold values v t h, with r a = 1.5 
andl/a„ = 1.8. Filled symbols correspond to the simulations, 
while open symbols are from an experiment on planar crack 
propagation in Plexiglas plates |6|, with (v) = 11 iim/s and 
a pixel size of 1.7pm 2 . The two data collapses (simulations 
and experiments) have been shifted on top of each other to 
facilitate comparison. 



cut-off in the distribution that scales as 



P(a) 



(6) 



with r a = 1.53 ± 0.05 and l/a v = 1.8 ± 0.1, see Fig. 
[2] There, we compare the simulation results with exper- 
imental data on planar crack propagation in Plexiglas 
plates @. Both sets of data can be collapsed by using 
exponent values r a 1.5 and l/a v « 1.8. By shifting 
the two data collapses on top of each other, we observe 
perfect agreement between the simulations and the ex- 
periment. Notice again that in Refs. jflH], the cut-offs 
of the distributions were included in the power law fits, 
thus overestimating the true exponent value. 



of the cluster size distributions scales with k as k 1 / <Tfc , 
with I/07. = 0.7 ± 0.08 for the AM approach, while we 
obtain l/a k = 0.8 ± 0.08 for the WTM. The different 
l/(7fc exponents have the same value within errorbars, 
and we suspect that the small variation observed might 
be related to a finite size effect. 



B. Comparison with experimental results 

In order to compare directly our results to the exper- 
iments @, we study WTM clusters as a function of a 
finite velocity threshold v t h- The threshold introduces a 



C. Link between global avalanches and local 
clusters 

To explain these results, we consider the break-up pro- 
cess of an avalanche into a number of clusters. When an 
avalanche grows by one unit (i.e. when a single element of 
the line moves forward by one unit), there are three pos- 
sibilities how this can affect the number N(s) of distinct 
clusters within an avalanche: i) a new cluster is initiated 
(N N + 1), ii) two existing clusters are merged into 
a single cluster (N — > N — 1), or iii) an existing cluster 
grows by one unit (N — > N). We associate probabilities 
pi, P2 and 1— p\— P2, respectively, to these events. These 
probabilities are dynamic variables as they depend on the 
instantaneous structure of the avalanche. On the average, 
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however, p\ = pi as p\ > pi would clearly indicate an in- 
crease in p 2 and a subsequent tendency for pi to decrease 
again: When the number of clusters grows, there is less 
and less space to create more clusters and the probabil- 
ity to merge existing clusters should increase. Similarly, 
pi < P2 would imply that sooner or later most of the 
clusters would have merged, and consequently the prob- 
ability P2 for further cluster merging decreases. Thus, for 
large enough avalanches, the number N > 1 of clusters 
during the growth of an avalanche follows a random walk 
up to s steps, with a reflective boundary at N = 1 (as 
each avalanche consists of at least one cluster) . Assuming 
for simplicity that the process is not correlated indicates 
that for s ^> 1, (N(s)) ~ s a , with a = 1/2. Our numeri- 
cal results confirm this relation, and a scaling form 



(N(s)) = s a N(s/k- 1/ak ) 



(7) 



with N(x) a scaling function, a = 0.47 ± 0.05 and a — 
0.76 ±0.08, see Fig. [5] (i.e. the cut-off scales with fc again 
with the same exponent l/cr^ ~ 0.75). Thus, the cluster 



sizes scale as a ~ s/(N(s)}, 
then implies that P{a) 



a 2 . P(a)da = P(s)ds 



with 



2r- 1. 



(8) 



In the present case, the value of the avalanche exponent 
r ss 1.25 yields r a w 1.5 for the clusters, in good agree- 
ment with our numerical and recent experimental results 
(i. 1!) . Notice that Eq. (J7J implies that s ~ a 2 holds only 
for s smaller than the cut-off scale. For larger avalanches 
one has a ~ s (as evidenced by the cross-over to a plateau 
in Fig. [5]), consistent with the observation that the cut- 
offs of the avalanche and cluster size distributions scale 
with the same exponent 1/crfe. Notice also that while de- 
riving Eq. (jSJ) we did not make assumptions about the 
precise form of the long-range kernel nor the spatial di- 
mension of the system, and we thus expect Eq. © to 
be valid for a wide class of avalanching systems where 
the avalanche dynamics is governed by long-range inter- 
actions. Notice that if the avalanche break-up process is 
correlated such that (N(s)) ~ s a with a ^ 1/2, Eq. © 
can be generalized to read r a = (r — a) /{I — a). 



D. Cluster morphology and roughness of the crack 

front 

In experiments @ and simulations Q it was found that 
the aspect ratio of the clusters was scaling as l y ~ , with 
£ ~ 0.6 and where l x and l v are the linear sizes of the 
cluster in the x and y directions, respectively. This re- 
sult is in agreement with the roughness of the crack front 
measured in early experiments, but disagrees with the 
known value of the roughness exponent of the crack line 
model ( = 0.39 [HI, l21| . Moreover, by investigating a 
larger range of length scales, a recent experimental work 
on roughness of the crack front has shown a crossover 
from a roughness exponent of around 0.6 to a lower value 
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FIG. 5: (color online) The scaling of the average number of 
clusters within a global avalanche of size s, for various values 
of k, and To = 1. The solid line is a power law fit of the form 
(N(s)} ~ s a , with q = 0.47. The inset shows a data collapse 
according to Eq. Q, with a = 0.47 and l/o k = 0.76. 



corresponding to the line model prediction of around 0.39 
[HI EH • This is surprising, since we have just shown that 
the experimental avalanche statistics is in perfect agree- 
ment with the predictions of the model. This puzzle can 
be resolved, noticing that Q = 0.39 is expected to hold 
only for length scales larger than the Larkin length L c 
[22| . which for our model scales like L c ~ Tq^/R 2 , where 
£ and R are the correlation length and the amplitude of 
the disorder (here £ = R = 1), respectively [l7[. Below 
L c , the crack line roughness should scale with an expo- 
nent Cl = 1/2. 

In the top panel of Fig. [6] we report the scaling of 
l y with l x for different values of Tq. By rescaling l x 
by L c ~ Tg, we observe a crossover for the scaling ex- 
ponent going from £l = 0.55 ± 0.05 at small scales to 
C = 0.39 ± 0.03 at large scales. In the top panel of Fig. [5] 
we also report the scaling of the root mean square height 
fluctuations (Ah(S) 2 ) 1/2 (with Ah(S) = h i+s - hi) of 
the crack front, expected to scale like (A/i(<5) 2 ) 1 / 2 ~ S^. 
Again, by rescaling 5 by Tq, we observe a crossover from 
a small scale regime with ^ = 0.48 ± 0.05 to the large 
scale value £ = 0.37±0.03. The bottom panel of Fig. [S] 
shows that similar conclusion can be made by consider- 
ing the power spectrum S(q) of the line profile, expected 
to scale as S(q) ~ q 2l ^ +1 . Rescaling the data by L c ~ Tq 
according to S(q) = T^S^Qq), we observe a crossover for 
the scaling exponent going from 2£l + 1 = 1-96 ± 0.05 for 
large q to 2£ + 1 rj 1.76 ± 0.04 for small q, in agreement 
with the analysis of root mean square height fluctuations. 
We have verified that in the present model, where the 
crack front is constrained to move along the y direction, 
there is no multiscaling p3j . Notice that while we ob- 
serve different crossover scales for the roughness of the 
line and the scaling of the morphology of the clusters, 
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FIG. 6: (color online) Top: The scaling of the aspect ratio of 
the local clusters (open symbols) and the root mean square 
height fluctuations (Ah(8) 2 } 1 ^ 2 of the line (filled symbols), 
for k — 0.0125 and different values of To. For large clusters 
their aspect ratio scales as l y ~ Zj, with £ = 0.39 ± 0.03 
(solid line), while smaller clusters are characterized by £l ps 
0.55 (dashed line). Also the roughness of the line, scaling 
like (Ah(S) 2 ) 1 ^ 2 ~ 5^, exhibits two distinct scaling regimes, 
with ~ 0.48 for small scales and £ = 0.37 ± 0.03 for large 
scales. Both sets of data have been collapsed by rescaling l x 
and S by the Larkin length L c ~ Tq. Bottom: The power 
spectra of the line profiles for different values of To. The main 
figure shows a collapse of the spectra, according to S(q) = 
rj)S'(rog), while the inset displays the unsealed power spectra. 
The power spectra scale as S(q) ~ g~' 2l = +1 ', with a cross-over 
separating regimes with 2C,l + 1 ~ 1.96 for short length scales 
and 2( + 1 ~ 1.76 for long length scales. 



both sets of data scale with the Larkin length L c , and 
thus the cross-over length is proportional to L c in both 
cases. However, the fact that these two length scales are 
not the same implies that there is a range of scales for 
which the roughness of the line appears to scale with the 
asymptotic exponent £ w 0.39, while the cluster aspect 



ratio is still characterized by the larger £l 0.55. It is 
probably this range of scales that has been investigated 
in earlier numerical studies 0]. 



IV. CONCLUSIONS 

To summarize, we have demonstrated how avalanches 
are broken into spatially disconnected clusters in systems 
with long range interactions. Both the avalanches as a 
whole and their localized parts or clusters exhibit scaling, 
but with different power law exponents, related by a scal- 
ing relation derived from a simple random walk argument 
for the avalanche break-up process. We also showed that 
large enough clusters exhibit aspect ratio scaling with an 
exponent consistent with the roughness exponent of the 
crack line model (13 . l2l| . Due to the general nature of 
both the model considered as well as our arguments, we 
expect these results to be applicable in a large class of 
avalanching systems with long range interactions, rang- 
ing from crack propagation to contact lines in wetting to 
domain walls in ferromagnets. 

Of particular interest here is our explanation of the dif- 
ferent scaling exponents observed for global avalanches 
and local clusters in the crack line model [7| , and its po- 
tential implications for further studies of avalanche phe- 
nomena in systems with long-range interactions. An in- 
teresting prospect would be to explore the possibilities 
to explain the observed acoustic emission exponents in 
peeling of paper sheets Q with similar arguments. The 
observation that the avalanche break-up process in the 
line model can be described by a simple random walk will 
be interesting to test in experiments on crack front prop- 
agation in Plexiglas plates, but also in other avalanching 
systems with long range interactions, such as for domain 
wall dynamics in ferromagnetic thin films [9j, as well as 
for imbibition [24| . 

An important point regarding the roughness of the 
crack front and the related cluster morphology is that 
their scaling properties depend on the range of length 
scales considered: The asymptotic scaling is observed 
only for length scales larger than a cross-over scale pro- 
portional to the Larkin length of the crack line. Earlier 
simulations 0] as well as experiments @ have most likely 
been probing the regime below this cross-over scale, pos- 
sibly explaining the different exponent values observed 
in those cases. However, for these small length scales 
additional complications arise due to the overhangs ob- 
served in the experimental crack profiles, which are ex- 
cluded by construction in the present line model. Such 
overhangs might be responsible for the multiscaling ob- 
served for the experimental crack profiles for short length 
scales [l6| • We are working to extend the crack front line 
model to include the possibility to form overhangs, and 
are planning to check if such an extension is sufficient to 
reproduce the experimentally observed short length scale 
scaling properties of the crack fronts. 

The connection between avalanches and clusters in 
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avalanching systems with long-range interactions opens 
up also interesting possibilities to understand the ob- 
served spatio-temporal correlations between avalanches 
in a number of systems, such as the space-time coupling 
of avalanches observed in plastically deforming crystals 
pEj , and even the fore- and aftershock sequences of earth- 
quakes: One possibility worth further studies would be to 
check if such correlated avalanches could in fact be local 
clusters, with the observed correlations arising naturally 



from the fact that the localized clusters are parts of the 
same global avalanche. 
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